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基于 /1 范 数 的 形状 快速 匹配 算法 
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(江南 大 学 物 联网 工程 学 院 , 江苏 无 锡 214122) 


摘 E: 针对 内 距离 形状 上 下 文 IDSC (inner-distance shape context) 和 轮廓 点 分 布 直 方 图 CPDH (contours points 
distribution histogram) 在 形状 相似 性 度量 中 直方 图 匹配 耗 时 长 ， 工 程 应 用 性 不 佳 的 问题 ， 提 出 了 一 种 用 EMD-Z1 测量 
轮廓 特征 直方 图 的 距离 的 方法 。EMD-Li1 在 原始 EMD (earth mover’s distance) 的 基础 上 融合 了 ZL1 范 数 ， 通 过 替换 地 面 距 
离 计 算 方法 ， 减 少 了 目标 函数 的 变量 ， 加 快 了 直方 图 匹配 的 速度 ， 能 够 快速 实现 形状 匹配 并 保持 较 好 的 检索 性 能 。 对 
形状 数据 集 进行 仿真 实验 的 结果 证 明 ， 该 方法 能 够 有 效 地 进行 数据 集 的 形状 识别 和 检索 ， 并 且 在 MNIST 数据 集 下 的 
匹配 速度 优 于 其 他 算法 。 
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Fast shape retrieval based on Li norm 


Wang Jianghui, Wu Xiaojun 
(School of IoT Engineering, Jiangnan University, Wuxi Jiangsu 214122, China) 


Abstract; In order to solve the problem that the histogram matching time is long and the engineering application is poor. This 
paper proposed to use EMD-L1 to measure the distance between two feature histograms. EMD-Li fusion of the Lı norm and the 
original Earth Mover's Distance (EMD) and replace the calculation of the ground distance. The algorithm reduce the number of 
unknown variables, could achieve shape matching quickly and has a good retrieval performance. With a great deal of experiments 
in several shape databases, the results show that the performance of novel method are superior to original algorithm. And the 
matching speed is better than other algorithms under the MNIST data set. 
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结构 6]。 基 于 轮廓 的 方法 有 多 尺度 法 ， 如 多 尺度 三 人 
面积 (triangle-area representation, TAR) i6, £ REHM 

形状 是 物体 的 一 项 重要 特征 ， 相 比 纹 理 、 灰 度 等 外 观 特 性 (multiscale convexity, MCC) [1 、 多 尺度 积分 不 变 
征 , 形状 具有 更 深层 次 的 视觉 信息 0 。 形 状 的 匹配 应 用 于 计 [8 等 。 这 类 方法 通过 构建 多 个 尺度 下 的 形状 特征 来 拓 
算 机 视觉 的 许多 领域 ， 如 目标 识别 ,图 像 检索 等 都 可 以 利 特征 信息 。 形 状 上 下 文 [9 (shape context, SC) 是 基于 轮 
用 形状 作为 分 类 与 识别 的 一 项 依据 。 其 中 形状 检索 又 可 以 分 廓 点 相对 位 置 关 系 的 经 典 方法 ， 该 算法 注重 轮廓 中 某 个 
为 2D 形状 与 3D 形状 检索 ， 因 为 最 近 几 年 深度 学 习 发 展 迅 点 与 其 他 所 有 点 的 空间 分 布 关 系 ， 但 是 其 方法 抗 噪 能 力 
速 得 到 了 众多 学 者 的 青睐 。 提 出 了 大 量 基 于 深度 学 习 的 3D 差 ， 采 用 薄板 样 条 函数 (TPS) 计算 直方 图 距离 耗 时 长 。 
形状 检索 算法 B4。 但 是 2D 形状 检索 依旧 是 计算 机 视觉 中 的 Ling 在 SC 的 基础 上 用 轮廓 点 的 内 部 距离 代替 欧 氏 距离 来 
代表 性 问题 。 二 维 形状 匹配 的 过 程 主 要 可 以 分 为 两 部 分 : 首先 衡量 轮廓 点 的 差异 ， 提 出 了 内 距离 形状 上 下 文 Cinner-distance 
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Ds 


设计 一 种 形状 描述 子 ， 即 通过 某 种 方式 来 描述 和 表示 形状 ; 然 shape context, IDSC) 09 方 法。 该 方法 对 非 刚性 物体 轮廓 的 形状 
后 计算 不 同形 状 的 描述 子 间 的 距离 来 衡量 形状 间 的 差异 完成 萌 述 具有 良好 的 表示 ， 并 采用 动态 规划 算法 度量 目标 形状 之 间 
匹配 。 该 方法 主要 分 为 基于 形状 组 成 点 集 、 形 状 边界 轮廓 的 距离 。 但 是 其 计算 复杂 度 高 ， 实 用 性 不 佳 。 

和 形状 内 部 骨架 3 类 来 描述 形状 特征 。 轮 廊 相 比 点 集 增加 Shu 和 Wut 提 出 了 一 种 2 维 形状 描述 符 一 一 轮廓 点 分 布 
了 点 的 顺序 信息 ， 骨 架 则 包含 了 形状 的 几何 特征 和 拓扑 直方 图 (contours points distribution histogram, CPDH). CPDH 
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以 目标 形状 的 最 小 外 接 圆 为 特征 提取 区 域 ， 以 轮廓 特征 点 分 布 。 心 圆 的 半径 , 8 为 区 域 Ri; 的 一 条 扇形 边 与 水 平方 向 ( 顺 时 针 旋 转 ) 
为 统计 对 象 ， 以 外 接 圆 圆 心 为 原点 建立 极 坐标 下 的 轮廓 点 分 布 。 所 成 夹 角 ，ni 为 落 在 区 域 Ri 中 的 点 的 个 数 。 这 样 一 幅 图 像 1 就 
直方 图 ,CPDH 用 地 球 移动 距离 (Earth Mover’s Distance, EMD ) 可 以 描述 成 R 个 三 元 组 所 构成 的 直方 图 有 厂 。 CPDH 具体 描述 细 
度量 直方 图 之 间 的 距离 。 但 是 CPDH 侧重 于 轮廓 局 部 特征 描述 节 参 考 文献 中. 

导致 其 在 形变 较 大 的 形状 数据 集 下 性 能 不 佳 ， 检 索 不 稳定 且 
EMD 算法 时 间 复 杂 度 高 ， 在 形状 匹配 时 速度 过 
本 文 提出 用 EMD-Li 度量 [9 由 轮廓 描述 子 生成 的 两 幅 形状 
轮廓 特征 直方 图 的 相似 度 ，EMD-LI 在 原始 EMD 包 基础 上 利用 
etn De dal E 
减少 函数 的 变量 和 约束 条 件 ， 降 低 了 计算 复杂 度 ， 加 快 
i oe 。 除 
此 之 外 工 : 范 数 在 特征 提取 中 也 有 去 噪 ， 减 少 异 常 值 的 优点 ， 可 
使 形状 描述 子 更 具 鲁 棒 性 通过 在 各 数据 集 下 检索 效果 的 实验 分 

析 ， 该 改进 算法 有 效 地 保持 了 IDSC 和 CPDH 在 良好 的 检索 性 1.2 IDSC 形状 描述 符 
能 ， 提 高 了 匹配 速度 ， 增 强 了 工程 应 用 型 。 首先 利用 二 阶梯 度 Canny 算 子 提取 目标 形状 的 部 分 轮廓 点 
= ee 集 Q={g,9;,4;,…,4,) ,qe IR*,n 表示 采样 点 数量 。 对 于 任意 采 
-i 样 点 都 可 以 用 其 余 n1 个 点 和 该 点 构造 的 向 量 表示 。 因 此 整个 
轮廓 特征 是 图 像 最 本 质 的 特征 ，CPDH 是 基于 轮廓 点 分 布 目标 形状 就 可 以 用 n*n-1 个 向 量 表示 。 这 些 向 量 包含 了 轮廓 点 
的 基础 上 提出 的 一 种 目标 特征 描述 符 。IDSC 是 基于 采样 点 之 在 整个 形状 中 的 距离 信息 和 方向 信息 。 方 向 参数 表示 把 一 个 
间 的 空间 分 布 关系 且 在 内 距离 的 基础 上 提出 的 一 种 具有 良好 表 ” 周 等 分 份 数 ， 距 离 参 数 是 等 分 半径 数 。 然 后 用 形状 统计 直方 
示 的 形状 描述 符 。 表示 得 到 的 向 量 集 。 
1.1 构建 CPDH 描述 符 内 距离 定义 为 形状 A 内 部 连接 两 点 m,n 的 最 短路 径 长 度 ， 
首先 采用 标准 的 Canny 算 子 对 图 像 目 标 轮廓 进行 提取 ,但 ] 4d (m,n; 人 A) 来 表示 。 最 短路 径 一 般 存 在 两 种 情况 : 若 点 m,n 
Ne Oo Ro eee 之 间 的 线段 完全 位 于 目标 形状 之 间 , 则 内 距离 等 价 于 欧 氏 距离 ; 
,采用 等 间隔 距离 方法 对 轮廓 点 进行 采样 ,得 到 轮廓 的 部 分 。 否则 应 用 Floyed-Warshall 最 短路 径 算 法 来 计算 每 一 对 采样 点 之 
. eee ek A E 间 的 距离 。 结 合 上 述 部 分 ， 可 得 到 目标 形状 的 IDSC。 如 图 3 所 
AMANITA IA 1 所 示 。 图 中 ，(a) 为 初始 图 像 ;(o) 提 取 的 mA, A 3.(a) 表 示 提 取 的 蝙蝠 形状 轮廓 并 标记 了 4 个 轮廓 点 ， 图 
采样 点 轮廓 图 像 ; (d) 为 目标 形状 的 最 小 外 接 圆 图 像 。 3.(b) 分 别 表示 相应 的 内 距离 形状 上 下 文 特征 。 
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图 2 极 坐标 下 轮廓 点 分 布 


E 


PR] 


— 


En 


(a) 原 图 像 人 b) 轮 廓 图 像 


(a) 形状 轮廓 标记 点 


IDSC at pt 1, bottom IDSC at pt 25, right 
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IDSC at pt 75, left 


(c) 采 样 点 图 像 (BSE IDSC at pt 50, top 
图 1 轮廓 提取 效果 图 
以 轮廓 点 集 P 为 目标 区 域 建立 其 最 小 外 接 圆 C, 以 圆心 (Cx， 


ey] 


Gy) 为 原点 建立 极 坐 标 系统 ,在 该 最 小 外 接 圆 内 以 同心 圆 方式 建 
立 如 图 2 所 示 扇 形 网 格 模型 来 统计 轮廓 点 的 分 布 特征 。 对 每 一 @) 标 记 点 对 应 的 直方 图 


个 扇形 区 域 Ri， 得 到 三 元 组 瓦 ={0,.0.n,} 其 中 : pi 为 第 i 个 同 图 3 形状 轮廓 和 统计 直方 图 
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录用 稿 


2 EMD 算法 


相似 性 度量 是 直接 影响 检索 效果 的 重要 环节 ， 人 们 一 
力 于 研究 如 何 用 数值 来 有 效 地 表示 图 像 在 特征 上 的 相似 程 
Rubner 等 人 0 提出 的 EMD 算法 广泛 应 用 于 线性 规划 中 的 


直 致 
度 。 


运输 


问题 。 直 觉 地 ， 本 文 将 两 幅 进 行 比较 的 图 像 看 做 是 图 像 特 
空域 的 两 个 不 同 分 布 ， 一 个 看 做 是 空间 散布 的 土 堆 ， 另 
做 是 土 坑 , EMD 算法 就 是 找 一 个 最 佳 路 径 来 完成 土 堆 的 运 


征 在 


一 个 看 


输 问 


题 ， 在 该 最 佳 路 径 下 所 做 的 功 可 看 做 是 两 个 分 布 间 的 距离 。 该 


算法 可 详细 地 形式 化 如 下 的 线性 规划 问题 : 
设 空间 SORE m NEH P =Ò Wa) nW) FB 
pi 个 土 扒 的 质量 为 ww ,同时 分 布 着 nm 个 土 坑 
Q =Ò WwW) (gw ,第 qi 个 土 坑 能 容纳 土 的 质量 为 Ww，， 
D=[dj] 表示 地 面 距离 EE 阵 。 需 要 找到 最 佳 运输 方案 下 =[ 方 ] 
这 里 fi eM HE pi 运输 到 土 坑 qi 的 土 的 质量 。 运 土 的 最 小 功 
可 定义 为 
Work( P,Q, F) = > di fy (1) 
i=l j=l 
该 式 需 满足 以 下 四 个 限制 条 件 : 
fi>0 1<i<ml<j<n (2) 
Sf <w, 1<i<m (3) 
j=l 
vf, sw, l<j<n (4) 
DÈS -miaf Sm | (5) 
i=l j=l i=l jal 
条 件 1 限定 只 能 把 土 从 土 堆 集 P 运 倒 土 坑 集 Q, 反 之 不 然 ; 
条 件 2 限定 土 堆 pi 最 多 所 能 运输 的 图 的 质量 不 超过 w,, ; 条 件 
3 限定 土 坑 q 最 多 所 能 承载 的 土地 质量 不 超过 wy Ph 4 表示 
尽 可 能 多 的 运输 土 。 只 要 该 运输 问题 得 到 解决 ， 就 找到 了 最 佳 


运输 方案 FE，EMD 距离 可 采用 最 终归 一 化 的 形式 : 
PIA 
EMD(P,Q)= =; (6) 
Lah, 
3 EMD-L1 度量 
CPDH 采用 原始 EMD 度量 轮廓 特征 直方 图 相似 度 ， 但 它 


存在 计算 复杂 度 高 , 对 微小 
卡 方 检验 计算 形状 直方 图 的 相似 度 得 到 代价 和 矩阵， 再 用 动 
JERE 


划算 法 (dynamicprogramming algorithm，DP) 得 出 代价 
匹配 。 虽 然 上 述 两 种 度量 算法 匹配 精度 高 但 是 计算 时 间 
不 利于 工程 应 用 。 针 对 这 一 不 足 ， 本 文采 用 EMD-L1 作为 


BAR SS PEPE BOSS GR RA © IDSC R) 
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中 最 
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EN. EMD-Li 算法 核心 在 于 用 Li 范 数 (曼哈顿 距离 》 葵 
原始 EMD 的 地 面 距离 。Li 范 数 〈 曼 哈 顿 距离 ) 是 指向 量 


代 了 
中 各 


个 元 素 的 绝对 值 之 和 ， 因 此 新 定义 的 地 


而 距离 如 式 (7) 所 示 。 


通 


过 利用 Li 范 数 取代 原始 地 
的 变量 个 数 , 加 快 了 直方 图 的 匹配 速 


面 距离 的 计算 方式 降低 了 


标 函 数 
于 Li 范 数 对 噪声 ， 


度 , 


异常 值 的 鲁 棒 性 优 于 Le 范 数 所 以 选取 Li 范 数 。 


di xı = 上 一 台 上 7 一 1 (7) 
其 中 : G ;deD 分 别 是 目标 形状 对 应 轮廓 特征 直方 图 下 的 块 。 
给 定 两 个 由 CPDH 生成 的 m 行 a 列 二 维 直方 图 P，Q， 分 
别 有 N 个 块 ，N =mxn， 作 如 下 的 定义 : 
a) 定 义 直方 图 块 的 索引 集 : 
7={G, filsis<mls jg<n},(ij) 表 示 直 方 图 中 的 一 个 块 
b) 定 义 块 之 间 流 的 集合 : v={(i,j,k,D):(i,j) ev,(k,D ev} 


c)F a) 可 得 ， 
Q={g;:(i,)) en} 


直方 图 P，Q 的 定义 : P={py:(i,j)en} , 


根据 上 述 三 个 定义 ， 重 新 定义 直方 图 


P 和 Q 间 的 EMD : 


OO) a PN DR hj (8) 
Fina 20 v(i, j,k 1) ev 
s.t. De = dqu V(k,1) en (9) 
ae fiixa=Py Vij)en 
F 表示 直方 图 P 到 Q 的 流 ; Suis ROR Gj) 到 (lol) 的 
流 ， 地 面 距离 dijxw 用 式 (7) 计 算 即 可 。 
为 了 进一步 优化 目标 函数 ， 将 v 划分 成 三 个 子 集 


= 每 一 
v= UV, Uv, 


个 子 集 对 应 一 和 


中 情形 下 的 流 : 


位 置 块 的 流 , 命 


vli jij) (Ljev) 表示 在 直方 图 同 


名 为 s-flows 

v ={(i, kD): (i jk, I) €v,d, jx1 
个 块 的 流 ,命名 为 nflows 
v ={(i j,k,1): (L441) E€ v,d; p41 > 


下 的 流 ,命名 为 fflows 
根据 以 上 定义 ， 可 以 将 人 flow & 


1) 表示 直方 图 中 相 邻 两 


1} 是 除 上 述 两 种 情况 番 


1 分 成 若干 个 n-flow。 这 是 


因为 在 整数 直方 图 块 上 Li 距离 可 以 


成 一 个 最 短路 径 系 统 。 例 


如 ， 给 定 一 个 fflow fijxni<k,j<l ,Li 地 面 距离 可 以 有 如 
下 分 解 : 
di jp di iia + ilikli T = 2 a, d, , x 十 2 dy; y+l,l (10) 
jsx<l 


综 上 所 述 ， 目 标 形状 P 与 查询 


EIR Q 对 应 CPDH 下 的 轮 


廓 直方 图 EMD-LI 距离 为 : 
DO it 
Eiji 29 V(i, j,k) €v; 
S.t. : 
Db! :(i, j,k.) € Vi (8i jaa = gi ) =b, v(i j)en 


(12) 
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录用 稿 Zi, $: 基于 工 1 范 数 的 形状 快速 匹配 算法 
其 中 : b; = Py— dy 表示 两 个 直方 图 在 块 (i,j) 的 轮廓 点 数 差 异 。 表 1 Kimia-25 数据 库 在 不 同方 法 下 的 检索 结果 比较 

标 函 数 11 和 约束 条 件 12 与 原始 EMD 相 比 ， 变 量 数 与 约束 方法 lst 2nd 3rd 
条 件 减少 。 另外 根据 Li 距离 的 特性 也 不 用 再 计算 块 与 块 之 间 的 SC 器 23 21 20 
也 面 距离 ， 很 大 程度 降低 了 算法 的 计算 复杂 度 。 CPDH+EMDI 25 24 22 

Sharvit et. all 23 21 20 
4 ”算法 时 间 复 杂 度 分 析 IDSC+DPI 25 24 25 
sie 0 ee a CPDH+EMD-L! 25 22 20 

改进 算法 中 的 优化 问题 将 目标 函数 的 变量 数 由 O(N?) ake BI) DSCEMDTI 35: 34 -34 
ON), N 表示 外 接 圆 划分 的 块 数 。EMD-L1 利用 L1 范 数 计 算 地 
面 距离 ， 地 提 于 [ 方 图 匹配 速度 。 原 始 的 EMD 算法 
nity oonan onia ayy CESTEESERTE 
ence a MAR ON), 通 过 实验 验 we 
证 改进 算法 的 形状 匹配 速度 得 到 了 有 效 提高 。 hw we mM 和 在 哥 可 妆 计 下 本 哥 刘 站 省 

NA 

5 验 结果 与 分 析 i | SS seSse22%4 64 
实验 结果 与 a 站 一 = - = 
mM a Toe 

KEWU wean nen dwn 

为 验证 本 文 所 提 方 法 的 有 效 性 , 在 不 同 数据 集 下 进行 实验 7 

p 、 、 、 (a)Kimia-25 形状 数据 集 (b)Kimia-99 形状 数据 集 
分 析 。 实 验 过 程 中 ， 分 别 将 数据 集中 每 一 样本 视 为 查询 形状 ， ee te 

oo ei peas JA ”两 组 形状 数据 集 剪影 
余下 样本 视 为 目标 形状 ,将 得 到 的 形状 距离 结果 按照 升序 排列 。 

如 果 匹 配 结果 显示 查询 形状 与 目标 形状 属于 相同 类 别 , 则 匹配 表 2 Kimia-99 数据 库 在 不 同方 法 下 的 检索 结果 比较 

结果 正确 ;否则 ,匹配 错误 。 方法 Ist 2nd 3rd 4th Sth 6th 7th 8th 9th 10th 全 部 
本 本 mee SC™ 97 91 88 85 84 77 75 66 56 37 756 

5.1 Kimie 25 = Kimla 22 MAARRE CPDH+EMD 96 94 94 87 88 82 80 70 65 22 808 

Kimia-25 数据 库 04 中 包含 6 个 类 别 共 25 个 样本 ， 样 本 示 IDSC+EMD-L1 99 99 98 98 96 90 85 82 76 67 893 
例如 图 4(a) 所 示 。Kimia-99 数据 库 05 包 含 9 个 类 别 各 11 个 样 CPDH+EMD-Li 99 91 85 85 79 80 73 78 68 53 791 
本 ， 样 本 实例 如 图 4.(b) 所 示 。 5.2 Articulated 形状 数据 集 

在 本 文 的 所 有 实验 中 ,为 确保 与 其 他 方法 比较 时 的 公平 性 ， Articulated 数据 集 00 包 含 8 个 类 ， 总 共 40 张 图 片 。 每 个 
设置 各 形状 轮廓 的 采样 点 数目 均 为 N=100。 在 Kimia-25 数据 ” 类 的 5 张 图 片 都 有 不 同 程度 的 形变 ， 如 图 6 所 示 。 该 数据 集 被 


A 
库 和 Kimia-99 数据 库 的 测试 过 程 中 , 目标 形状 最 小 外 接 圆 圆 周 FA PURE, EIB AE BARRA. E 
12 等 分 ， 圆 心 角 半径 5 等 分 ， 即 Spx126= 60 个 区 域 块 。 也 可 以 检验 算法 对 形变 的 鲁 棒 性 。 该 实验 采用 交叉 验证 ， 分 别 
在 Kimia-25 数据 库 进 行 仿真 实验 的 结果 见 表 1, 表 1 列 出 用 卡 方 检验 x o EMD-L 度量 同 一 形状 描述 子 生成 的 直方 图 之 
了 不 同方 法 在 Kimia-25 数据 集 下 的 检索 结果 , 给 出 了 与 查询 形 | 间 的 相似 度 。 将 返回 与 查询 形状 最 相似 的 1~4 组 目标 形状 且 分 
状 最 相似 的 1~3 组 目标 形状 中 分 类 正确 的 具体 数目 。 从 表 1 的” 类 正确 的 具体 数目 作为 评价 算法 检索 性 能 的 标准 。 分 析 表 3 的 
结果 中 可 以 看 出 ,本 文 所 改进 的 IDSC 与 CPDH 效果 略 优 于 经 实验 数据 , 第 一 组 与 第 二 组 实验 结果 相 比 , 说 明 EMD-Li H 7? 
方法 SC, 主要 原因 在 于 SC 在 构建 每 个 轮廓 点 的 局 部 直方 图 ”在 形变 较 明显 的 数据 集 下 的 检索 效果 更 好 ， 和 鲁 棒 性 更 加 ; 分析 
时 ， 若 半径 过 短 ， 则 局 部 直方 图 不 具备 良好 的 鉴别 能 力 ; 若 半 ”第 二 组 与 第 三 组 的 实验 结果 ,IDSC 对 形状 的 鉴别 能 力 优 于 SC, 
径 过 长 ， 则 局 部 直方 图 会 包含 较 多 的 噪声 。 但 是 所 提 算 法 检索 分析 第 一 组 和 第 五 组 实验 ,可 以 看 出 CPDH 对 形状 形变 的 描述 
效果 弱 于 CPDH+EMD, 与 IDSC+DP。 经 分 析 LI 距离 蔡 代 原始 ”能 力 更 优 于 SC 因为 IDSC 对 形状 的 描述 中 包含 了 轮廓 点 的 相 
EMD 地 面 距离 时 忽略 了 轮廓 点 之 间 的 空间 信息 , 导致 检索 效果 ”对 位 置 关 系 且 SC 算法 抗 噪 性 能 差 。 
略 有 下 降 。 5.3 MPEG-7 数据 集 
在 Kimia-99 数据 库 进 行 仿真 实验 的 结果 见 表 2， 表 2 分 别 根据 本 文 提出 的 算法 ， 选 用 MPEG-709 数 据 集 进行 仿真 实 
给 出 了 与 查询 形状 最 相似 的 1~10 组 目标 形状 分 类 正确 的 具体 验 ， 该 数据 集 共 有 1400 个 形状 剪影 ， 均 分 为 70 个 类 别 ， 每 个 
数目 及 合计 结果 。 从 表 中 可 以 发 现 ,本 文 算法 检索 性 能 优 于 SC, 类 别 中 共有 20 个 样本 。 
改进 的 CPDH 算法 弱 于 CPDH+EMD。 但 是 检索 到 数据 集 第 8 ”5.3.1 CPDH+EMD 与 CPDH+EMD-L1 的 比较 
类 至 第 10 类 时 , 本 文 算 法 比 另 外 两 种 算法 更 加 稳定 , 返回 有 效 K 4 284  CPDH+EMD, CPDH+EMD-L; 算法 的 比较 ， 本 
的 形状 数目 更 多 。 原因 在 于 采取 LI 距离 对 微小 形变 的 鲁 棒 性 更 “次 实验 在 MPEG-7 Set B 的 子 图 像 库 中 进行 , 该 子 库 包括 216 幅 
好 。 另 外 改进 后 的 IDSC 算法 检索 性 能 明显 优 于 SC, 形状 图 像 ， 分 为 18 个 类 别 ， 每 类 12 幅 。 从 表 4 的 实验 结果 可 
CPDH+EMD. .证 明 EMD-L1 能 很 好 的 保持 IDSC 对 形状 的 描述 以 看 出 ， 虽 然 CPDH+EMD 返回 的 有 效 形 状 总 数 多 于 
能 力 。 CPDH+EMD-LI， 但 是 分 析 第 8 类 到 第 11 类 的 检索 结果 ， 尤 其 
是 第 10 类 至 第 11 类 EMD 返回 的 形状 数量 下 降 较 快 而 EMD- 
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ER, $: 基于 Ll 范 数 的 形状 快速 匹配 算法 


Li 则 平稳 减少 且 第 11 类 目 远 多 于 原始 算法 ， 说 明 
F EMD-Li 的 CPDH 在 类 数 较 多 的 情况 下 对 形变 的 鲁 棒 性 更 
但 是 由 于 Li 范 数 的 特性 忽略 了 轮廓 点 之 间 的 空间 信息 导致 总 


体 效 果 不 如 原 算法 。 
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图 5 Articulated 形状 数据 集 


表 3 Articulated 形状 数据 集 在 不 同方 法 下 的 检索 结果 比较 


方法 Topl Top2 Top3 Top4 

SC+EMD-L1 36/40 19/40 10/40 8/40 

SC+X20 20/40 10/40 11/40 5/40 

IDSCT+X2 40/40 34/40 35/40 27/40 

IDSC+EMD-L1 39/40 3940 34/40 32/40 

CPDH+EMD-Li 40/40 12/40 14/40 10/40 

表 4 MPEG-7 子 库 中 的 检索 结果 比较 

方法 Ist 2nd 3rd 4th 5th 6th 7th 8th 9th 10th 11th 总 数 
CPDH+EMD"! 214 215 209 204 200 193 187 180 168 146 114 2030 
CPDH+EMD-Li 216 206 199 185 186 178 176 167 154 144 135 1946 


每 个 待 检索 形状 都 要 与 数据 集中 所 有 形状 比较 ， 取 实 


验 结果 中 与 每 个 待 检索 形状 最 相似 的 40 个 形状 ， 统 计 这 40 个 


形状 中 与 待 检索 形状 同类 别 的 数目 ， 并 计算 对 应 的 


查 全 率 。 从 


表 5 可 以 看 出 EMD-L1 提升 了 IDSC 的 检索 性 
数 形状 检索 算法 。 证 明了 本 文 算 法 的 有 效 性 。 
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表 5 各 算法 检索 性 能 比较 

方法 检索 精度 (Bullseye) 
CSSD3 75.44% 
Visual Parts!'7! 76.45% 
SC+TPS(y’)"! 76.51% 
Curve Edit!" 78.17% 
Distance Set!”°! 78.38% 
Gen.ModelP20 80.03% 
IDSC(y’) 85.40% 
IDSC+EMD-L1 86.56% 


6 ”旋转 不 变性 测试 


超过 了 大 多 


为 了 验证 IDSC+EMD-L1 的 旋转 不 变性 。 本 次 实验 从 


MPEG-7 数据 集中 随机 选取 10 
每 一 幅 图 像 幅 分 别 沿 着 顺 时 针 和 逆 时 针 方 向 依次 旋 


30° ,45° , 60° ,75” 得 到 12 幅 新 的 图 像 。 包 括 10 幅 基 准 图 


画图 像 作 为 基准 图 像 ， 然后 将 这 


转 5° 15°, 


象 ， 整 个 旋转 数据 集 有 130 AK. 
图 像 作 为 查询 图 像 的 检索 结果 ， 若 有 检索 错误 的 图 


5.3.2 IDSC+EMD-L1 与 其 他 方法 的 比较 


背景 颜色 设置 为 灰色 。 从 表 6 可 以 看 出 检索 返回 结 


表 6 列 出 了 将 10 基准 幅 
像 将 其 表格 


果 只 有 两 幅 


图 像 检 索 错 误 ， 其 余 118 图 像 检 索 正确 ， 说 明 改 进 


本 实验 采用 MPEG-7 数据 集 并 使 用 Bull-eye07 方 法 确定 检 然 保持 良好 的 旋转 不 变性 。 
表 6 IDSC+EMD-L1 旋转 不 变性 测试 结果 
查询 图 像 4 5 8 
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后 的 算法 仍 
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速度 测试 


本 实验 采用 两 组 算法 进行 对 比 ,第 一 组 实验 : CPDH+EMD- 
L1, CPDH+EMD, SC; 第 二 组 实验 : IDSC+DP，IDSC+EMD- 
Ll1，SC。 采 用 的 数据 集 是 MNIST[231 手 写 数 字 图 像 库 的 训练 图 
BE, 共有 10000 WAR, 包含 0 到 9 共 10 种 数字 。 从 中 随机 
选取 1000 幅 进行 实验 (数字 形状 复杂 度 从 低 到 高 ， 由 此 计算 出 
的 识别 速度 更 具有 一 般 性 )。 实 验 结果 如 图 6、7 所 示 。 实 验证 
明 改 进 后 算法 的 运行 时 间 显著 优 于 原 算法 与 SC, 更 适合 于 工程 


应 用 。 
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图 6 第 一 组 算法 时 间 测 试 
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7 第 二 组 算法 时 间 测 试 


8 RE 


本 文 基于 CPDH Fil IDSC 分 别 用 EMD-L; #4% EMD 和 DP 
算法 作为 轮廓 点 直方 图 的 度量 准则 ， 其 主要 优点 在 于 Li 范 数 
(曼哈顿 距离 》 对 于 形变 有 较 好 的 鲁 棒 性 ， 而 且 通 过 减少 目标 

函数 的 变量 数 降低 了 算法 的 计算 复杂 度 。 从 实验 结果 来 看 ， 所 
进 算 法 有 良好 的 检索 能 力 ， 加 快 了 形状 匹配 速度 。 由 于 最 近 
几 年 深度 学 习 大 热 ， 以 神经 网 络 为 基础 做 3D 形状 检索 将 是 作 
者 下 一 步 研究 的 方向 。 
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